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ABSTRACT 

Warped accretion disks have attracted intensive attention because of their critical role on shaping the spin of 
supermassive massive black holes (SMBHs) through the Bardeen-Petterson effect, a general relativistic effect 
that leads to final alignments or anti-alignments between black holes and warped accretion disks. We study 
such alignment processes by explicitly taking into account the finite sizes of accretion disks and the episodic 
lifetimes of AGNs that delineate the duration of gas fueling onto accretion disks. We employ an approximate 
global model to simulate the evolution of accretion disks, allowing to determine the gravitomagnetic torque 
that drives the alignments in a quite simple way. We then track down the evolutionary paths for mass and 
spin of black holes both in a single activity episode and over a series of episodes. Given with randomly 
and isotropically oriented gas fueling over episodes, we calculate the spin evolution with different episodic 
lifetimes and find that it is quite sensitive to the lifetimes. We therefore propose that spin distribution of 
SMBHs can place constraints on the episodic lifetimes of AGNs and vice versa. Applications of our results on 
the observed spin distributions of SMBHs and the observed episodic lifetimes of AGNs are discussed, although 
both the measurements at present are yet ambiguous to draw a firm conclusion. Our prescription can be easily 
incorporated into semi-analytic models for black hole growth and spin evolution. 

Subject headings: accretion, accretion disks — black hole physics — galaxies: active — quasars: general 


1. INTRODUCTION 

Accretion disks surrounding supermassive black holes 
(SMBHs) are believed to be the source engines of enormous 
radiative p ower of active galactic nuclei (AG Ns) in the uni¬ 
verse (e.g.. lSalpeterlll964tlLvnden-Belllll969l) . The shapes of 
accretion disks—flat or warped—depend on the orientation of 
the angular momentum of gas fueling from the circumnuclear 
region with respect to the rotating orientation of the central 
black hole. A wide variety of evidence have indicated that in 
many circumstances these two orientations are misaligned. 


First, the orientation of AGNs seems uncorrelated with 
that of their host galaxies, well established by observa- 


tions 

Kinney et al.l 2000; 

Schmitt et al.: 2003; Gallimore et al. 

2006 

Munoz Marin et al. 

2007; Shen et al.l 201T&|Lagos et al. 

2011 

and also confirmed by numerical simulations (e.g.. 

Hopkins et al.l 12012 and references therein). This suggests 


that gas channeled onto the central SMBH either loses 
the memory of its initial direction in the h ost g alaxy or 
comes from outside t he galaxy (e g.. iDubois et al.1120121 
iNavakshin et al.ll2 012t iGaspari et a ll (20131) . most likely due 
to galaxy mergers lKendall et al.ll2003h . Secondly, the current 
hierarchical framework of galaxy formation and evolution 
predicts that repeat galaxy mergers, in addition to secu lar evo¬ 
lution, trigger SMBH activities (e.g jBenson & Boweili201()l) 
and gas accretion plausibly proceeds at episodic and random 
phases (in particular for minor mergers; e.g.. lKing & Pringlel 
120061: I Wang et aLll2006l) . Thirdly, recent studies on SMBH 
spin through quantifying the radiative efficiency of AGN pop¬ 
ulations under the thin accretion disk model showed that 
SMBH spin undergoes a decline w ith cosmic time sine z ~ 2 
(IWang et all 120091 : iLi et all[201 2l) . potentially implying that 
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(partially) random accretion takes place (iKing et all 12008 ). 
Such evolution ary b ehavior has also been reproduced by 
IVolonteri et all d2013l) in their semi-analytic model of SMBH 
growth and evolution (see their Figure 10). 

Once gas fueling is inclined with respect to the central spin¬ 
ning black hole, a warped accretion disk forms inevitably. The 
Lense-Thirring torque arising from the frame dragging effect 
will lead to precessions of the inclined accretion disk and the 
black hole around each other. Detailed dynamics of warped 
accretion disks has been studied extensively through theo- 
retical analys i s (e.g., iPringlel 1 1992k iPapaloizou & Linl [ 1995b 
lQgilvie| |1999t iLubow et all 20021) a nd numerical simulations 
(e.g. Nelson & Papaloizou 1999L I2000t iLodato & Pringlel 


120061120071 iFragila120091 IZhuravlev et aUl2014l) . The inten- 

sive applications of warped accretion disks lie at their critical 
role on shaping the spin of the central SMBHs thro ugh the 
Bardeen-Petterson effect dBardeen & Pettersonll 19751 ). In this 
effect, the presence of strong viscosity, a most like case for 
normal AGNs, damps out the differential Lense-Thirring pre¬ 
cession and induce the inner portion of the disk to align or 
anti-align with the spin axis of the hole. Accompanied with 
warp propagation, the entire disk will tend to align or anti¬ 
align with th e hole dScheuer & Feilerlll996t IKing et aUl2005l : 
IMartin et~akll2007l) . 

Previous studies on the alignments between black holes 
and their warped accretion disks mainly focus on two as¬ 
pects. The first assumes that accretion disks have infinite 
extensions and constant profiles. Black holes are always 
driven toward t he disks and the b l ack hole-disk systems end 
up alignments (iScheuer & Feilerl 119961 IMartin et alJ 120071: 
iPerego et al.ll2009h . On the contrary, the second assumes that 
accretion disks are finite in both size and mass and the total 
angul ar momentum of the systems is conserved dKing et all 
12005112008; ). which means that there is no gas supply to the 
disks. Such systems can end up a configuration of either 
alignment or anti-alignment, depending on the initial inch- 
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nations and the ratios of the angular momenta between disks 
and black holes dKing et al.ll2005l) . 

In reality, it is well known that accretion disks are trun¬ 
cated by their own self-gravity in the outer regions where it 
becomes dominated over the gravity from black holes (e.g., 

Paczvnskill978llShlosman & Begelmanll987l : ICollin & Zahnl 

19991) . Outside the truncation radius, the disks fragment 
into dense clum ps and finally turn to be a hotbed of star 
formations (e.g.. iGoodmanl 12003b iLevinl 120071 IWang et ahl 
120101 LTiang & Goodmanll201 lh . It comes naturally that the 
stellar radiation and supernova explosions feedback to the 
ambient gas, exciting turbulence that acts to transfer the 
angular momentum and leads to gas supply to the accre- 
tion d i sks from the further oute r regions d Wada & Normanl 
20021 iKawakatu & Wadal 120081 iKumar & Johnsonl 120101 
Wang et al.ll2010l) .Therefore, :o appropriately study warped 
accretion disks, it is very practical to take into account both 
the finite sizes limited by the self-gravity and the plausible 
gas feedin g at the outer boundaries (see also the discussions 
of lLi et ak 120131) . The lifetime that the gas feeding lasts for, 
in principle, can be determined from the observed episodic 
lifetime of AGNs, for which there are several techniques, al¬ 
though the results are still primitive fsee iMartliiil120041 for a 
review). 

In this paper, we revisit the alignment processes of black 
holes with their warped accretion disks by including the fi¬ 
nite sizes of accretion disks and the continuous gas feeding 
controlled by the episodic lifetimes of AGNs. For the conve¬ 
nience of numerical implementation, we introduce a param¬ 
eter “the lifetime of gas fueling” to describe the duration of 
gas supply. We employ an approximate global model fol¬ 
lowing jKumaretai] (120081) to describe the accretion process. 
This enables us to study the alignments of finite-size accre¬ 
tion disks in a quite simply way, avoiding being immersed in 
the time-consuming numerical simulations on the evolution of 
accretion disks. 

The paper is organized as follows. Section 2 describes the 
properties of warped accretion disks and Section 3 derives 
the basic evolution equations for the black hole-accretion 
disk system. In Section 4, we present the results for 
alignments/anti-alignments and spin evolution in a single ac¬ 
cretion event and over a series of events. We then show 
the implications of our calculations for spin distribution and 
episodic lifetimes of AGNs in Section 5. Discussions on plau¬ 
sible caveats of our calculations and conclusions are given in 
Section 6. 

2. PROPERTIES OF WARPED ACCRETION DISKS 

In the weak-field limit, the frame-dragging effect due to a 
spinning black hole exerts a Lense-Thirring _toK]ue on each 
annulus of the accretion disk dBardeen & Pettersonll975l) . By 
denoting the angular momentum of the annulus at a radius R 
as L(R ) and that of the black hole as ./|,. the Lense-Thirring 
torque reads 

, 2 GJ h xL(R) 

^=-2 (1) 

where G is the gravitational constant and c is the speed of 
light. This Lense-Thirring torque leads the annuli to pre- 
cess around the rotating axis of the black hole at a rate of 
2GJh/c 2 R } . Since such differential precessing rate rapidly 
declines with radius as R 3 , the inner region of the accre¬ 
tion disk is intensely affected in the beginning and gradu¬ 
ally warps propagate outwards over the entire disk. Detailed 


theoretical analysis showed that the dynamics of warped 
accretion disks is characterized into two regimes depend¬ 
ing o n the importance of the pressure force and viscous 
force (iPa paloizou & Lin! 11995 ). If adopting the standard a- 
prescription of viscosity, warps propagate as bending waves 
for inviscid or sufficiently thick disks (a < H //?), whereas 
propagate in a diffusive fashion for strongly viscous disks 
(a > H/R). Here H is the semi-thickness of disks. For sim¬ 
plicity, we only consider thin and viscous accretion disks with 
Keplerian rotation, in which a > H/R is generally satisfied. 
This is the most likely case for accretion disks in AGNs with 
the Eddington ratio at around 0.01-1, the largest population 
in AG N surveys (e.g., iMcLure & Dunlopl l2004t IShen et ahl 

mu. 

In the presence of high viscosity, the differential precess 
in the inner disk will be damped out rapidly, giving rising 
to an inner fla_t_region (so called Bardeen-Petterson effect; 
iBardeen & Pettersonlll975 T This region extends to the ra¬ 
dius where the precessing timescale is comparable with the 
warp-diffusion timescale. The precessing time scale reads 
c 2 R 3 /2GJh and the warp-diffusion time scale reads R 1 /vo, 
leading to the warp radius as 

2GJ h 

R w-y, (2) 

vic- 

where v 2 is the vertical shear viscosity governing the warp 
diffusion through the accretion disks. 

Integrating over all the annuli yields the total torque exerted 
on the accretion disk 


4t tG 


./h X 


'L(R) 

R 2 


dR. 


(3) 


As a reaction, the black hole, in the meantime, suffers an equal 
but opposite torque that cause its angular momentum to pro¬ 
cess as well. By taking into account the angular momentum 
carried by mass accretion, the angular momentum of the black 
hole evolves following 

^ = jf h U h -r L x, (4) 

dt 

where ;V/ h is the mass accretion rate onto the black hole and 
£ ms is the specific angular momentum at the marginal stable 
orbit R ms . Hereafter a hat symbol on top of a vector represents 
the corresponding direction vector. 

Similarly, the angular momentum of disks changes with 
time as 

^ = -MhfmJh + T LT +jf, (5) 

dt 

where J<\ is angular momentum summed over all the annuli 
and Jf is angular momentum change due to gas fueling. Note 
that 7 lt is always orthogonal to ,/ h , therefore, it can be ex¬ 
pressed in a form of (iKing et al.ll2005l) 


Tut = K\Jh x ,/ d + /G./h x (,/h x J,i), (6) 


where the first term in the right hand side corresponds to the 
precession and the second term corresponds to change in the 
angle between the hole and the disk. We keep in mind that 
both K\ and /G are generally p ositi ve as shown in our previous 
numerical analysis iLi et al.ll2013l) . 

Considering that the inner region of the disk is aligned with 
the black hole and the Lense-Thirring torque drops off rapidly 
with radius, the major contribution to T\-\ in Equation (0 
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comes from the integral around the warp radius R w . We there¬ 
fore make an approximate for as 


K 2 


2t tG J h L(R w ) 
R 


(7) 


Below we will show that such an approximation of Kj yields 
results that exact ly agree with the previous analytical studie s 
for small warps (iScheuer & Feilerl 1 1996b IMartin et al.ll2007l) . 
Note that K\ only governs the precession of the black hole- 
disk system and therefore is not important for the present pur¬ 
pose. 


3. BASIC EVOLUTION EQUATIONS 
3.1. Accretion Disks and the Inclination Angles 

From above equations, we can now determine the evolution 
of the angular momenta Jh and / d and the angle 6 between 
the black hole and the accretion disk. We first define a life¬ 
time of gas fueling A7f as during which the accretion disk is 
continuously fed at a mass rate limited by the Eddington rate 
Me dd = L Edi /ric 2 = 2.2M & r/Z\M Q yr _1 , where L Edd is the Ed¬ 
dington luminosity, = Mh/ 1 O 8 M 0 , and r/ = 0.1 r/_i is the 
radiative efficiency for mass accretion (see below Section [L2l 
for definition). 

Note again the inner region of the disk is aligned with the 
black hole. Multiplying Jh to Equation 0 yields 

-j- = M d £ ms , (8) 

dt 

and multiplying J d to Equation (0 yields 


where 9q is the initial value of 9{t) at t = 0 and T = Jh/Ki = 
c 2 R w /2nGL(R w ). 

For small warps, we extend the trigonometric functions in 
Equation (fl2l > around 9 ~ 0 and it is trivial to show that the 
corresponding solution decays exponentially as 

6{t) = 6 > 0 exp ■ (13) 

This is the analytic solution fou nd inlS cheuer & Feilerl (119961 
see their Equation (11)) 0 and IMartin et all (120071 see their 
Equation (52)) and the time scale of the exponential decay T 
exactly agrees with their derivations. 


3.1.2. Phase II: Gas Fueling Quenched 

After gas fueling is quenched, i.e., J\ = 0, the total angular 
momentum of the system (black hole and accretion disk) be¬ 
comes conserved. The disk subjects to both alignment and 
viscous diffusion. The former decreases the angular mo¬ 
mentum of the disk (/Q > 0 in Equation 0), whereas the 
later transfers the angular momentum outwards. Substituting 
Equations 0 and 0 into Equation ( ITOl ) and with some sim¬ 
ple mathematical implementations, we arrive at 


d cos 9 
dt 


M\\ ^rns 

Jd 


sin~ 9 + —— sin” ( 

Jh 


1 + —cos 9 
Jd 


(14) 


Following lKumar et all (120081) . we employ an approximate 
global model to simulate the influence of mass accretion on 
the disk properties so that accretion rate M can be determined 
in a quite simple fashion. To this end, we define the effective 
radius of the disk R d afl 


—- = -M h l ms cos9-Kisin 2 9+JfCosd, (9) 

dt 


Jd_ 

M d 


■■£(R d ) = (GM h R d ) 1 / 2 . 


(15) 


where we assume that the angular momentum of gas fueling 
is aligned with J d . The angle between the black hole and the 
accretion disk changes with time as 


dcos9 dj h ■ J d dJ d dJ h 

JdJh — 7 — =- 7 -y h cos 6 »—-/ d cos 6 >—. ( 10 ) 

dt dt dt dt 


We distinguish the evolution of 9 into two phases as below. 


3.1.1. Phase I: Gas Fueling 

During the phase with gas fueling, the accretion disk is con¬ 
tinuously fed so that it maintains a preferred steady angular 
momentum distribution. As a result, one shall expect that the 
term due to gas fueling cancels out the other terms in Equa¬ 
tions 0 and 0 and J d and its magnitude J d remain somewhat 
constant, i.e., dj d /dt = 0 and dJ d /dt = 0. With Equations (3) 
and (7), we can simplify Equation (ITOl) into 


d cos 9 
dt 



( 11 ) 


Since K 2 is positive, 9 always tends to decrease towards 
zero regardless of its initial value, i.e., alignmenLbetween 
the disk and the b l ack hole (see alsolScheuer & FeileJ 1 19961 

IMartin et al. 1120071 iPerego etal. 1120091k 

In this phase, as the disk maintains its profile steadily, K 2 
can be deemed to be a constant. There exists an analytic solu¬ 
tion for Equation (fTTb 


1 - cos 9 1 - cos 9 d / 2t \ 

1 + cos 9 1 + cos 9q 6 ^ \ T ) ’ 


(12) 


The mean global accretion rate is then given by 


M h = - 


dM d M d 

dt Urn 


where the accretion time scale reads 


vdRiY 


(16) 


(17) 


where u\ is the normal viscosity associated to the azimuthal 
shear due to the differential rotation of the disk. With the 
standard a-prescription,;/] is expressed by 

iq =afl K H 2 = ah 2 n K R 2 , (18) 


where li = H/R is the disk’s aspect ratio and U K is the Keple- 
rian angular momentum. Combining above equations, we can 
obtain 


M h = 


ah 2 G 2 M 2 M 4 d 

4 


(19) 


This equation apparently indicates the influences of alignment 
and accretion on M d . From Equation 0, alignment processes 
reduce the disk angular momentum / d at a rate of ~ K 2 but 
leave the disk mass M d unchanged, therefore, increasing the 


3 Note that iScheuer & Feilej (1996ft omitted a factor of \/2 in their Equa¬ 
tion (11). 

4 We stress that the definition of aims to simulate the disk evolution in 
a global way. It has the same order of magnitude as, but is not necessarily 
identical to, the disk’s physical outer truncated radius. 
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Fig. 1.— Evolution of inclination angle 0 between black holes and disks for different lifetimes of gas fueling A7f. The initial values are (Left) 6 = 30° and 
(Right) 6 = 150°, respectively. Open circles in left panel represent the analytic solution for small warps: an exponentially decay with an e-folding timescale 
T = 2.2 x 10 5 yr (see Equation CCD). 



Time (10 6 yr) 



Fig. 2. — Same as Figure[T|but for m ass a ccretion rate. The initial mass accretion rate is M\j! = 0) = 0.25MEdd- In left panel, open circles represent the 
asymptotic evolution following Equation ED In right panel, mass accretion rates are limited by the Eddington rate. 


accretion rate ;V/ h . On the contrary, mass accretion reduces 
the disk mass but approximately keeps the disk angular mo¬ 
mentum unchanged, considering that angular momentum loss 
due to the accretion inflow through the marginal stable orbit 
is usually negligible, therefore decreasing M h. 

Indeed, if only considering the influence of accretion rate, 
there exists an asymptotic solution to Equation © (see 
iKumar et al.ll2008h 


M d (f) = M d (f = 0) 


1 + 


3 1 

— 

*acc J 


-1/3 


and 


M h (t) = 


M d (t = 0) 


/' 

acc 


i + 


3 1 

— 

*'acc J 


-4/3 


( 20 ) 


( 21 ) 


where f' cc = f acc (0) is the accretion time scale at t = 0. For 
small warp cases, we expect that the mass accretion rate ap¬ 
proaches the above asymptotic behaviors, i.e., M d oc r 4 3 for 

t > 4c- 

With mass accretion rate, we now rewrite Ki in a more con¬ 
venient form. Previous studies showed that the vertical vis¬ 
cosity v 2 is related to the normal azimuthal viscosity u\ by 


(lOgilviell 1 999b iLodato & Pringlel [20071 

120091) 

Vl _ fv 

v\ 2a 2 ’ 


see also lPerego et al.l 

( 22 ) 


where /„ is a coefficient of order of unity. As a zeroth-order 
approximation, we neglect the detailed shape of warped ac¬ 
cretion disks and make use of the estimate for surface den¬ 
sity E r; Mh/37rj/i and for the angular momentum L(R) r; 
M d sJGMhR as in flat disks. This is reasonable since 
accretion disks reach maximally warped around the warp ra- 
dius R w and then rapidly approach flat geometry beyond R w 

dScheuer & FeilerlFl996fc lMartin et al.l[2007l : iChen et al.1120091 

iLi et a I. II 201 3h. Combining Equations (O and Q gives 


K 2 » ^ VGM h R w . (23) 

For small warp amplitude, f v = 1, whereas for large warp 
amplitude, /„ moderat ely decreases below unity, dependin g 
on the w arp am plitude (lOgilviell 1999tlLodato & Pricell2010h . 
lOgilviel (119991) developed a prescription to determine /„ for 
any warp amplitude. However, this needs to know the de¬ 
tailed warp profile as a function of radius, therefore beyond 





































LI et al. 


5 




Time (10 6 yr) 


Fig. 3.— Same as Figure[T]but for mass growth of black holes. The initial mass of the black hole is Mh(t = 0) = 10 8 Mq. 




Fig. 4.— Same as FigureQ]but for spin evolution of black holes. The initial black hole spin is a = 0.5. 


the scope of the present work. We simply set /„ = 1 through¬ 
out our calculations and below we will show the reliability of 
our results in light of this simplification. 


3.2. Black Hole Evolution 


The direction of the inner flat region of the disk depends 
on the inclinations of the incoming gas flow and the black 
hole: if 9 < 7r/2, it aligns with the hole; while if 0 > 7t/ 2, 
it anti-aligns with the hole. By taking into account the mass 
loss carried away by radiation, one obtains that the black hole 
grows as 

^ = (1 -rf)M h , (24) 

at 

and by combining with Equation ©, the black hole spin 
evolves as 


da M h 

— = (1 -m — 
dt ' M h 


1 4 

1-tiRgC 


■ —2a 


(25) 


where r] is the radiative efficiency depending on the spin pa¬ 
rameter a via 

and the specific angular momentum is calculated as 
„ /- 1 =F2ar _3 / 2 + a 2 r“ 2 

£ ms = ± \ZGM h R ms ^_ 3i _- l±2ar - 3/2)1/2 - ( 2? ) 


for a < 1 and 

/-1 ± r -1 / 2 + r~ l =F r -3 / 2 

4is = ± VGM h R ms (i _|_2r _ i/ 2 )i/ 2 ’ 

for a = 1, where r = R ms /R g , R« = GM h /c 2 is the gravitational 
radius, and the upper sign refers to prograde accretion while 
the lower sign refers to retrograde accretion dBardeen et all 


4. RESULTS 

According to the standard accretion disk model, the disk’s 
aspect ratio h = H /R is in an order of 10 -3 ~ 1CT 2 and weakly 
depends on radius in proportional to R i22i> , on black hole 
mass in proportional to M h ^ 10 , and on mass accretion rate 

in proportional to M^/ 5 dSh akura & Sunvaevlll973l: s ee also 

Collin-Souffrin & Dumontill990l : iNataraian & Pringle! 1 19981 

King et al.ll2008h . An accretion disk becomes self-gravitating 
when its mass reaches a factor h of the black hole mass 

5 Note that, according to the calculation of lThornel H9741) , the radiation 
emitted by the geometrically-thin disk and swallowed by the hole produces a 
counteracting torque that prevents the hole span up beyond an upper limit of 
a ~ 0.998. In our calculations, the accretion may proceed from retrograde to 
prograde fashion over episodes. This effectively means than the black hole 
spin may suffer a flip from a = — 1 to a = 1. For the sake of simplicity, we 
neglect the counteracting torque and assume that black holes can carry any 
spin between -1 and 1. 
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(e.g.. lPringldll98l]) . We therefore adopt the initial disk mass 
Md ~ hM tl . The initial effective radius of the disk Rd is deter¬ 
mined by the required initial mass accretion rate and accretion 
time scale given with h. To be specific, from Equation (IT9l >. 
Mh can be alternatively expressed through Rd and h by 


M h 

^Edd 


= 0.25 (——) ( — 
V0.1/ V0.1 




5 x 10“ 3 

- 3/ V M h 8 1 


2.1xl0 2 R g J \10 *M q 
The corresponding accretion time scale is 

. -i / h 


face = 8.8 x 10 5 yr(^-) 


Rd 


3/2 


5 x 10 

( M h 


1-3 


2.7x10 2 R g J \10 S M Q J' 


(29) 


(30) 


This is the typical timescale for accretion disks truncated by 
self-gravitational collapse (see, e.g. JPerego et al .1120091) . The 
ratio between the angular momenta of the disk and the hole is 


Jd 

Jh 


0.26 a 


Rd 


5 x IQ” 3 / V 2.7 x 10 3 /?„ 


1/2 


(31) 


Our calculating procedure is as follow. We start with a 
black hole that has a mass of Mh and a spin of a. The accre¬ 
tion disk is assigned an initial mass of Md = /)M h and angular 
momentum of Jd = M^GM^R^) 1 ^ 2 . The initial inclination an¬ 
gle between the black hole hole and the disk is 9. We then 
evolve Equations ([9| and (fl6l> for the angular momentum and 
mass of the disk, Equations (l24l > and < 125 b for the mass and 
spin of the black hole, and finally Equations (1111 and (ITU for 
the inclination angle. The mass accretion rate Mh is deter¬ 
mined by Equation (IT9l >. We keep in mind that in the phase 
with gas fueling, the disk’s mass and angular momentum are 
unchanged. For a series of accretion events, we just repeat the 
above procedure but with the mass and spin of the black hole 
inherit from the previous episode. 

We assume that in the beginning of an episode, the disk has 
already been formed and neglect the stage for the formation 
of the disk. It is expected that during the formation stage, the 
accretion rate must be very insignificant so that the influence 
on black hole spin evolution is negligible. Also, we implicitly 
assumed that the time lag between the successive episodes is 
long enough to let the disk be consumed on the accretion time 
scale before the next episode starts. This can be verified by the 
fact that the observed “duty cycle” is generally m uch less than 
unity (e.g., IShankar et al.ll2009t iLi et al.lf2012l) . Duty cycle 
measures the fraction of active black holes to total black holes, 
and equivalently indicates the time fraction of active phases 
oyer total lif etime (active and inactive) for a single black hole 
dWang et al .1120061) . 

In following calculations, unless stated otherwise, we set 
a = 0.1, h = 5x 10 -3 , R d = 2.7 x 10 3 M~ 2/3 r] 2 _{ 3 R g , M h = 
10 8 Mq, and a = 0.5 as the fiducial values. It is worth men¬ 
tioning that the choice of Rd dependent on black hole mass 
and the radiative efficiency aims to let the dimensionless ac¬ 
cretion rate Mh/ Me dd at abo ut 0.25, the typical value observed 
in luminous AGN su rveys (iKollmeier et al.l120061 IShen et akl 
l2008t ILi et a 1.112001) . We also set an upper limit of mass ac¬ 
cretion rate by Mh < Me dd considering that beyond the Ed¬ 


dington limit, radiative feedback may be important to regulate 
the accretion rate at around the Eddington limit. 

4.1. Alignments and Anti-Alignments 

In Figure [I] we demonstrate evolution of the inclination 
angle between the disk and the black hole for different life¬ 
times of gas fueling: A7f = (0,1,2,4,10) x 10 5 yr. The ini¬ 
tial value of 9 is set to be 9q = n/6 and 9q = 5tt/6 in the 
left and right panels, respectively. As can be seen, for the 
case of 9q = 7t/ 6, the inclination angle rapidly decreases on 
a timescale of 10 5 yr, which is th e typical estimate i n previ¬ 
ous n umerous literature ( see, e.g. JPerego et al.l2009l:lLi et al.l 
120131 and references therein). As the lifetime of gas fueling 
A7f increases, the decay rate of 9 is slightly reduced. This is 
because during the gas fueling, the disk maintains its orien¬ 
tation and only the black hole is driven to align to the disk. 
While after the gas fueling is quenched, both the disk and 
the hole change their orientations towards the direction of the 
total angular momentum of the system. This can be further 
verified from Equations (fill and (ITTl i. Interestingly, the be¬ 
havior of 9 for A7f = 10 6 yr, namely, the disk is continuously 
fed throughout the time range of cal culations, is exactly iden¬ 
tica l to the analytic solu tion found in lScheuer & FeileTI (119961) 
and lMartin et al.l (120071) . 

The evolution of 9 is much more complicated for the case 
of 9q = 57t/ 6 as shown in the right panel of FigureQ] With no 
gas fueling (A7} = 0), the inclination angle rises up to 7 r all the 
way. Indeed, the instantaneous behavior of 9 is governed by 
the ratio Jd/Jh in Equation ( 1 1 4| i: for Jd/Jh > - cos 9 , the incli¬ 
nation 9 decreases, while for Jd/Jh < -cos 9, the inclination 
9 increases (see also lLi et al.li2013l) . As A7} increases, the fi¬ 
nal configuration of the system transits from anti-alignment 
(9 = 7 r) to co-alignment (0 = 0). We note there is a steep 
decreasing/increasing trend of 9 just after the gas fueling is 
quenched, ascribed to the rapid reduce of the angular momen¬ 
tum of the disk and hence enhancement of mass accretion rate 
due to the viscous dissipation arising from alignments. See 
below for detailed explanations. 

Figure[2]shows the corresponding mass accretion rates. For 
the sake of comparison, the asymptotic evolution of accre¬ 
tion rate by Equation (fTil is superposed in the left panel. If 
there is no gas fueling, as the inclination angle tends to zero, 
the angular momentum of the disk approach constant. As a 
result, the mass accretion rate evolves asymptotically follow¬ 
ing Equation (l2ll . In the right panel, the initial value of 9 
is 9q = 5n/6. Instead of a gradual decay, there exists a no¬ 
table peak of accretion rate, which even reaches the Eddington 
limit Me dd in some cases. The similar behavior of an enhance¬ 
ment of the mass accretion rate also appears in the previous 
numerical simulations o n warped accretion d isks (see, e.g., 
iLodato & Pringlell2006l and iNixon et al.l 120121) . This can be 
easily understood in terms of Equations ©, (1X31 ). and (IT6l) : 
strong dissipation arising from alignments leads to an intense 
reduce of the magnitude of ./ cl ; the effective radius Rd accord¬ 
ingly shrinks to maintain a Keplerian rotating disk. As a con¬ 
sequence, the accretion time scale f acc is shortened and mass 
accretion rate Mh is significantly enhanced. In reality, emer¬ 
gence of the enhancement of accretion rate depends on com¬ 
petition between the viscous-dissipation driven deduce of Jd 
(due to warp alignments) and mass-accretion driven deduce 
of Md- This is the reason why there is no accretion rate en¬ 
hancement for some curves in Figure [2] 

By assuming the angular momentum conservation of the 
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M h (10 B M q ) 

Fig. 5.— Evolution of spin magnitude with mass growth during a series of 
accretion events. The initial spin is a = 0.9, 0.5, and 0.0 for red, blue, and 
black lines, respectively. The gas fueling lifetime is A If = 10 6 yr for upper 
three lines, and there is no gas fueling for lower three lines. 

system, IKing et aTl ( 2005 ) derived a condition for occurrence 
of anti-alignments: the initial angle between the disk and 
the hole satisfies cos (A) < This equation is inten¬ 

sively used in subseq uent studies on cosmological evolution 


of SMBH spins (e.g., Volonteri et al. 20071 Kina et al. 

2008; 

Lagos et al.l 

2009; Fanidakis et al. 2011c Dotti et al. 

2013; 

Dubois et al. 

20141 Sesana et al.l 20141). We relax the assump- 


tion of the conservation of angular momentum and let the ac¬ 
cretion disk steadily fueled for a lifetime of A7]. Since the to¬ 
tal angular momentum of the s ystem is not conserved, it is in¬ 
appropriate to directly use the IKing et al.l d2005l) ’s condition. 
In the right panel of Figure [Q the initial ratios A/A = 0.40 
and cos^o = -0.87, therefore the condition cos 0q < ~A /27 h 
is well established. However, for A7f >2x 10 5 yr, the system 
ends up an alignment instead of an anti-alignment. 

Nevertheless, the ratio A/A remains an important indicator 
of alignments or anti-alignments (if there is no gas fueling). 
Small ratios mean that the angular momentum of black holes 
dominates over that of disks, therefore it is comparatively eas¬ 
ier to drive the disks towards alignments or anti-alignments. 
Conversely, if the angular momentum of the disks is domi¬ 
nated, the holes are “lighter” to change their orientations. This 
is somehow similar to the cases that the disks are continuously 
fueled: the ho les tend to align t ow ards the disks (see also th e 
discussions of lKing et al.1120051 and iLodato & Pringlejl2006l) . 

4.2. Spin Evolution 
4.2.1. A Single Accretion Event 

In Figures [3 and |4] we show evolution of the black hole’s 
properties, i.e., mass and spin, respectively, for initial inclina¬ 
tion angles 0o = 7t/ 6 and Oq = 5n/6 with a set of gas fueling 
lifetimes A7). Compared with the case of Oo = ~/6, we can 
find that the black hole mass growth for 0q = 57r/6 is much 
more efficient. This is because the Eddington limit is larger 
with lower radiative efficiency due to retrograde accretion for 
0 > 7t/ 2 and also because the accretion rate is enhanced as 
plotted in Figure [3 The spin evolution in Figure [4] clearly 
illustrates the important influences of gas fueling. As the 
lifetime of gas fueling increases, the black hole is gradually 
driven to align with the disk regardless of the initial inclina¬ 
tion. As a result, in right panel of Figure Q] there exists a 
trend for A 7) >2x 10 5 yr that spin first declines due to retro¬ 
grade accretion and then increases when the inclination angle 


transits to 0 < 7t/2. In addition, it can be found that the spin 
changes are relatively more significant for black holes with 
A7f < 2 x 10 5 yr, because these holes undergo retrograde ac¬ 
cretion throughout the episode, which carries larger angular 
momentum compared with the prograde accretion (see Equa¬ 
tions (EUi and @). 

4.2.2. A Series of Accretion Events 

The spin of black holes will reach an equilibrium value after 
a series of accretion events, de pending on the fraction of pro¬ 
grade and retrograde accretion (IKing et ali2008lfSesana et~ahl 
2014). Evidently, this fraction is intimately relevant to the de¬ 
grees of anisotropy in the orientations of accretion disks with 
respect to the holes. If all the accretion events occur with the 
initial inclination angles 0q <n/2, the black holes accrete in a 
coherent fashion and their spin will be spun up all the way to 
the maximum (a = 1). On the other hand, the inclination an¬ 
gles Op may distribute randomly oyer the whole space (0 ~ n) 
(e.g., IKing et al.ll2008k IWang et aPl2009t iHopkins et al.ll2012t 

iLi et al.ll2012l) . For simplicity, we only focus on the random 
distribution of Oo in proportional to sin6*o over (0 ~ n), so that 
both alignments and anti-alignments can occurs in our cal¬ 
culations. We first generate a random inclination angle and 
then apply the prescription described above to track down the 
mass growth and spin evolution of black holes. The calcula¬ 
tion ends when the accretion rate decreases to one tenth of its 
initial value, i.e., 0.025Me<m- 

In Figure |3 we illustrate how black hole spin evolves along 
their mass growth during a series of accretion events. We 
stress that we adopt the effective radius of accretion disks de¬ 
pendent on black hole mass as Rd = 2.7 x 10 3 A7 8 2 3 r/j| 3 R g in 
order to guarantee the initial dimensionless accretion rate at 
0.25 in Equation ( l29l >. Once there is no gas fueling, irrespec¬ 
tive of the initial spin values, black holes lose the memory of 
their initial status and approach a relatively low spin roughly 
after doubling their mass. Then as mass growth, there is a 
weak trend that their spin gradually falls off. Let’s give an in¬ 
tuitive explanation for such weakly negative dependence of 
spin on mass through the ratio A/A = L Since we adopt 

Rd oc , Equation (OTT) yields a ex A/ h 1 ' 3 . Of course, the 
real situation is much more complicated, but the negative de¬ 
pendence is retained (see also below for further discussions). 

Once the disks are fueled for a lifetime longer that the typ¬ 
ical alignment timescale (e.g., A7f = 10 6 yr), as expected, the 
holes are always aligned to the disks and are rapidly spun up 
to fast rotating. However, we note that the final spin never 
reaches the maximal value a - 1 because during the align¬ 
ments, there exists a period with retrograde accretion that 
spins down the holes. This indicates that even though the fi¬ 
nal configuration in an accretion event is alignment instead 
of anti-alignment, the spin decrease during the time with ret¬ 
rograde accretion still contributes to the final spin status and 
therefore should be appropriately included. 

5. IMPLICATIONS FOR SPIN DISTRIBUTION AND 
EPISODIC LIFETIMES OF AGNS 

5.1. Monte-Carlo Simulations 

We perform Monte-Carlo simulations to explore the spin 
distribution with the above described prescription. The free 
parameters a, h, and Rd are set by their fiducial values as in 
Section [4] Initially, we assign to a black hole a mass drawn 
from a uniform distribution in logarithm over (10 6 ~ 10 8 )M Q 
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log M h (M q ) 




log M h (M Q ) 


FlO. 6.— Spin distribution with black hole mass from Monte-Carlo simulations for different gas fueling lifetime A7f-: (a) 0, (b) 10 5 yr, (c) 10 5 M g 1 yr, and (d) 
10 5 A/g ,5 yr. Negative spin mean s that black holes a re u ndergoing retrogr ade accretion. Data poi nts w ith error b ars in panel (c) are the measured spin through 
broad Iron Ka line compiled bv IBrennemanl 120131) and [Reynolds! 12013 1 (see also Table 2 of lSesana et al.ll2(TT4| ). Horizontal dotted lines represents black hole 
spin a = 1 and -1, respectively. 


TABLE 1 

Episodic Lifetimes of AGNs Constrained from Observations in the Literature. 


Object Redshift Episodic Lifetime (Myr) Telescope Reference 


Proximity Effect 


Q 0302-003 

3.285 

10-30 

VLT 

1,2 

HE 2347-4342 

2.885 

~ 25 

VLT 

3 

KP 76 

2.466 

> 25 

Keck, HIRES 

4 

KP 77 

2.535 

16-33 

Keck, HIRES 

4 

HS 1700+6416 

2.748 

>20 

HST/COS 

5 

Quasar Pair Sample 

~ 2.2 

~ 1 

Keck, Lick, NOAO, SDSS 

6 



Fluorescent Ly a 

emission 


HS 1549+1919 

~ 2.84 

> 1-3 

Keck 

7 

QSO 0420-388 

~ 3.1 

>60 

VLT 

8 

Quasar Sample 

~ 2.7 

1-20 

Keck 

9 


REF ERE NCES. — (1) Worseck & Wisotzki <2006); (2) Jakobsen et al. (2003); (3) Worse ck et all 
(2007); (4) Goncalves et al. (2008); (5) McOuinn & Worseck (2014); (6) Kirkman & Tvtler (2008); (7) 

lAdelberger et alU2006l): (8)ICantalupo et al] <20071); (9) Trainor & Steidel 120131). 


and a spin drawn uniformly over (0 ~ 1). For this newly 
generated black hole, we first let it grow with 100 accre¬ 
tion episodes, to remove the memory of its initial condi¬ 
tions. Afterwards, this hole continues to accrete for a series 
of episodes. We randomly assign the number of episodes be¬ 
tween (0 ~ 900) so that in average the hole undergoes 500 
times activities. We then record the final mass and spin infor¬ 
mation of the hole. In each accretion episode, the inclination 
0o between the hole and its accretion disk follows a uniform 
distribution of cos 9q over (-1 ~ 1). An individual event is re¬ 


garded to be terminated once the accretion rate decreases be¬ 
low a tenth of its initial value. We repeat such procedure 500 
times and are finally left with a sample of 500 black holes. 

In Figure [6] we plot the spin distribution with black 
hole mass for different gas fueling lifetime A7) = 0, 10 5 yr, 
1() s A7 g 1 yr, and 10 5 Mg 5 yr in panels (a)-(d), respectively. In 
panel (a), there is no gas fueling and the disks will be con¬ 
sumed within an accretion timescale f acc ~ 8.8 x 10 5 yr (see 
Equation (l30l i). As the ratio Jd/Jh decreases with black hole 
mass, the accretion disks of massive black holes are easier 
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log M h (M 0 ) 

Fig. 7.— Radiative efficiency with black hole mass for gas fueling lifetime 
AT, = 10 5 Mg' 5 yr. Negative spin means that the black holes are undergoing 
retrograde accretion. Gray points w ith erro rbars are the radiative efficiency of 
local quasars estimated bv IDavis & Laorl 1201 1 ) and shaded area represents 
the radiative efficiency at z ~ 1 estimated by [Li et al. 1 2012). Horizontal 
dotted lines represents r] = 0.42, 0.057, and 0.038 for black hole spin a = 1, 
0, and -1, respectively. To guide the eye, a solid line with a slope of 0.5 is 
plotted in the bottom right corner. 

change their orientations and thus anti-alignments become 
more probable. This results in a negative correlation between 
spin and mass. A similar trend was also found by dDotti et al] 
I2Q13L corresponding to F = 0 in their model). When we switch 
on the gas fueling for a time of A7} = 10 5 yr in panel (b), black 
holes are always driven to align towards the disks during this 
period. Therefore, the fraction of retrograde accretion is sig¬ 
nificantly reduced and all the black holes carry spin \a\ > 0.5. 

We also choose a mass-dependent lifetime of gas fueling 
for the sake of comparison with observations. In panel (c), 
the gas fueling lifetime is A7f = 10 5 Mg'yr, indicating that the 
disks surrounding 10 7 M q black holes will be fueled for 10 6 yr, 
far longer than the alignment timescale. Accordingly, smaller 
black holes are overall spun up to maximum rotating. There 
is a strong anti-correlation between spin and mass. On the 
other hand, in panel (d), we reverse the mass dependence of 
fueling lifetime to A7f = 10 5 Mg 5 . As expected, massive black 
holes at ~ 10 9 M 0 have maximum spin (a = 1) whereas lower 
massive black holes in this case are not affected (since 7f is 
much less that the accretion timescale and can be neglected) 
and just follow the distribution in panel (a). 

5.2. Comparison with Observations 

We now compare our calculated spin distributions with 
observations. Notwithstanding measuring spin of massive 
black holes is yet challenging so far, there are sev eral de- 


veloping techniques shedding useful light (e 

g., 'Wang et al.! 

2009; Czerny et al. 201 1; Davis & Laor 2011 

Li et al. 2012; 

Brenneman 2013; Done et al. 2013tlRevnolds 

20lflWu et al. 

2013|). Among these techniques, we focus on the direct mea- 
surements using the broad Iron Ka line fitting (Brenneman 


120131 Reynolds! 120131 and references therein) and the indi¬ 


rect estimates using the radiative efficiency of thin accretion 
disks as a surrogate of the black ho le spin (see Equation (26l ) : 
iDavis & Laorll20lH iLi et al. 112012l) . Both these two methods 
had been applied to a sample or a population of AGNs, allow¬ 
ing us to perform straightforward comparison. Unfortunately, 
each method bears a large source of uncertainties that are not 
well understood and the obtained results seem incompatible 
at this stage. Therefore, we perform comparison separately as 
follows. 


In panel (c) of Figure[6l we superpose the spin of 19 sources 
measured from Iron K a lines, most of which are Seyfert I or 
narrow line Seyfert I AGNs res i ding in spi r als or lenticular 
galaxies (see lSesana et alJf2014l) . iRevnoldsI (120 13h proposed 
that in this sample low massive black holes (~ 10 7 Mq) are 
spinning fast while black holes with mass > 10 8 M Q may have 
moderate spin. Our calculations show that a lifetime of gas fu¬ 
eling as A7f = 10 5 Mg 1 can generally reproduce such observed 
spin distribution. This means low massive black holes require 
longer gas fueling to reach alignments and to accrete coher¬ 
ently. 

On the other hand, IDavis & Laorl ( 201 l b reconstructed the 
spectral energy distributions for a sample of 80 Palomar- 
Green quasars and estimated their individual radiative ef¬ 
ficiency of accretion flows based on the thin disk model. 
They found a tight correlation between the efficiency with 

1 /2 

black hole mass app roximately as 77 oc M h ' (but see also 
iRaimundo et al.l2012l) . plausibly i mplying a raise in the black 
hole sp in with mas s. Meanwhil e, ILi et al.1 (120121) applied the 
ISoItanl s argument (iSoItanll 19821) to solve the continuity equa¬ 
tion of SMBH demography in active and quiescent galaxies 
(ILi et all 120111 ) and quantified the radiative efficiency with 
redshift and black hole mass. Their results suggested that the 
efficiency undergoes cosmological evolution and generally 
there is an increasing trend o f the efficiency with black hole 
mass at redshift z > 1 (s ee also lCao & Lil2008l : IVolonteri et alJ 
120131 lUeda et all 120141) . Figure [7] compares the radiative ef¬ 
ficiency from our calculations by adopting the gas fueling 
lifetime A7} = 10 5 Mg 5 yr (as in panel (d) of Figur e [6jl with 
these from observational constraints of IDavis & Laorl (1201 ll) 
and lLi et~afl (120121 at z ~ 1). Here, the radiative efficiency is 
calculated from black hole spin as Equation (l26l > by assum¬ 
ing that the thin disk model applies. As stated above, in the 
context of chaotic accretion without gas fueling, high massive 
black holes (~ 10 9 Mq) have low spin because of net spin- 
down by retrograde accretion with anti-alignments. There¬ 
fore, to maintain the observed high radiative efficiency (and 
high spin), sufficient prolonged gas fueling is in need for these 
holes. 

Longer gas fueling means that there is plenty of gas reser¬ 
voir progressively channeled into the accretion disks and the 
black holes are accordingly long-lived. There are several tech¬ 
niques developed for measu ring or estimating episodic life¬ 
times of AGNs, see iMartinil (20041) for a review. In Table Q] 
we summarize the observational constraints on episodic life¬ 
times of AGNs through the proximity effect and the fluo¬ 
rescent Lyo emission. Note that both these constraints are 
presently only realizable for high-redshift, luminous quasars, 
indicating that there reside high massive black holes (plau¬ 
sibly > 10 9 Mq). It seems that the observed AGN episodic 
lifetimes are generally ~ 10 7 yr but with large uncertainties. 
Such long lifetimes strongly imply that high massive black 
holes are rapidly spinning at high redshift, consistent with 
previous studies probed through the radiative efficie ncy (e.g., 
IWang et all2009tlLi et al J2QI 2llTrakhtenbrotl2014l) . Our cal¬ 
culations in the present work further show that if the correla¬ 
tion between the efficiency and the black hole mass is reliable, 
low massive black hole must accrete in more chaotic way and 
a gas fueling lifetime as A7f ~ 1 0 5 Mg 5 can roughly reproduce 
the observed slope 77 oc M h ' . This gives rise to an episodic 
lifetime of > 3 x 10 6 yr for SMBHs with Mh > 1 0 9 Mq, coin¬ 
cident with the observations. 
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6. DISCUSSIONS AND CONCLUSIONS 

The main differenc es of the present work compared with the 
previous studies (e.g JPerego et al.12009!iDotti et al.l2013l) lie 
at: 1) we include the influences of continuous gas fueling on 
the alignments between black holes and their accretion disks; 
and Tf we consider the finite-size accretion disks (see also 
iKing et akll2008l) . There are several simplifications in our cal¬ 
culations meriting further explanations. First, we neglect the 
detailed structures and properties of accretion disks and em¬ 
ploy an approximate global model to simulate the evolution 
of disks under alignments and accretion (iKumar et al .112008b . 
We are thus unable to handle the angular momentum trans¬ 
portation throughout the disks. This may be important when 
the disks suffer large warps (9 ~ 7t/2) in which non-linear 
effects cause the disk structures around the warp radius to 
change dramatically and the simple global approximate plau¬ 
sibly becomes invalid. However, since the spin parameter in¬ 
vokes the whole accretion history, we expect that our results 
are insensitive to this effect. Secondly, the main free param¬ 
eters in our calculations are the disk’s aspect ratio h and the 
effective radius R^. In terms of the mass accretion rate, h 
and R ( \ is degenerated along lines of constant h 2 R^ 1 . As men¬ 
tioned above, the standard accretion disk model shows that 
h is almost insensitive to radius and only weakly depends on 
black hole mass and ma ss accretion rate dNataraian & Pringle! 
119981 IKing et al.ll2008l) . We therefore fix h to the intermedi¬ 
ate value between 10 -3 ~ 10 -2 . The effective radius is chosen 
by a specified mass accretion rate in Equation (l29l >. which is 
roughly at ~ 0.25MEdd from AGN surveys. We recalculate the 
spin distribution for h between 1 0 -2 and 10 -3 and find the re¬ 
sults are qualitatively unchanged. Thirdly, the Lense-Thirring 
torque K 2 calculated b y Equation (l23l> is also a well approx¬ 
imate for small warps dLodato & Pringlej|2006l) . When disks 
are strongly warped, on one hand, the relation M ~ 3ni /\£ no 
longer applies, in particular in the region around the warp ra¬ 
dius; on the other hand, the vertical shear viscosity V 2 deviates 
from the relation in Equation (l22l> and turns t o be decreasing 
as wa rp amplitude increases (Qgilvie 1 999r. iLodato & Pricel 
120101) . To see how our results depend on f v , we artificially 
set an extreme value of /„ =0.1 for large inclination angle 
9 £ (7r/3-27r/3) and show the comparison with /„ = 1 in Fig¬ 
ure [8] The ch oice of a such range is based on the study by 
iPerego et all ( 2009 ) who found that the non-linear term be¬ 
comes important for angular momentum transportation when 
the inclination angle > 7t/3. Black hole spin magnitude is sys¬ 
tematically smaller because of the reduced torque K 2 for large 
warps. However, we can find that the changes are moderate 
with such an extreme /„ value and the overall spin distribution 
with mass is preserved. 

For simplicity, we also assume that disk fueling is chaotic 
so that the initial inclination angle 9q between the disks and 
black holes distributes randomly over (0 ~ 7 r). iDotti et ahl 
(120131) proposed a comprehensive modeling on the anisotropy 
of 9q using a parameter F to describe the fraction of accre¬ 
tion events with 9q > n/2. Therefore, F = 0.5 stands for 
fully ch aotic fueling and F = 0 stands for fully coherent ac¬ 
cretion. ISesana et all (120141) further extended this prescription 
and connected the anisotropy parameter F with the dynamic 
properties of the host galaxies. They postulated that the ratio 
of the rotation velocity to the velocity dispersion of the galax¬ 
ies mirrors the degree of anisotropy of the gas component 
channeled into the nuclear regions. They argued that the rota¬ 
tion velocity measures the bulk rotation of the galaxies while 



Fig. 8.— Influence of fv on spin distribution. Blue points represent the 
case /„ = 1 for all the 6 s and red points represent the case fv = 0.1 for 8 £ 
(7t/ 3—27r/3). The gas fueling lifetime is ATf = 10 5 Mg- 5 yr. 


the velocity dispersion measures how chaotic the galaxi es ar e. 
This seems in line with the recent finding of lHo & Kim (120141) 
that the virial factor in reverberation mapping, which encodes 
the geometric and kinetic properties of the broad-line regions 
of AGNs, is related with the large-scale morphology of the 
host galaxy bulges. However, numerical simulations showed 
that the angular momentum of gas inflowing into the nuclear 
regions (^pc) will lose its mem ory of the initial direction 
on the the larger sc ales (~ kpc; iBarnes & Hernquistl 1 1 996b 
iHopkins et al.ll20f2h . This is well supported by the observa¬ 
tions that AGNs are misaligned with the host galaxies (e.g. , 
KinngY_gtalJ j2000l: iGallimore et all 120061 : IShen et all 120101 : 
Lagos et aD 120111) . These conflicts reflect our poor under¬ 
standing on the detailed processes in black hole fueling from 
the galaxy scale to the central disk scale. We await future ob¬ 
servations to justify these conflicts and we content with our 
present assumption. 

We consider only thin accretion disks with the disks’ aspect 
ratio II /R <C a to guarantee that the viscosi ty is str ong enough 
to validate the Bardeen-Petterson effect dPanaloizou & Linl 
1199511 . Therefore, our approach only applies to the accretion 
disks with dimensionless accretion rate at 10 -2 ~ 1, 

which are believed to correspond to the thin disk regime. Be¬ 
yond this range, accretion disks transit to be thick and the 
Bardeen-Petterson effect may be weak or even disappear, as 
confirmed by recent m agnetohydrodynamics (MHD) numer - 
ical simulations (e.g., iFragilel 120091 : IZhuravlev et all 120141). 
Howev er, the heuristic MHD simulations bv lMcKinnev et all 
( 2013 1) revealed a new “magneto-spin alignment” mechanism 
that works in thick disks with strong magnetic. Unlike the 
Bardeen-Petterson effect, this mechanism arises from the 
magnetic torque of the black hole’s magnetosphere, which is 
already aligned with the black hole spin axis due to the frame¬ 
dragging forces. It is unknown yet how the spin evolves for 
thick warped disks and more investigations are highly worth¬ 
while. 

In conclusion, we study the alignments of black holes and 
their accretion disks by taking into account the finite sizes of 
disks and the continuous gas fueling. Our results show that, 
with fiducial values for the free parameters, the lifetime of gas 
fueling is crucial to the alignments/anti-alignments and there¬ 
fore to the spin evolution. By applying our prescription to a 
series of accretion activities and assuming that the disk orien¬ 
tations are chaotically distributed over episodes, we compared 
our calculated spin distribution with the spin measurements 
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through Iron K a line fitting and through the radiative effi¬ 
ciency and made attempt to place constraints on the lifetime 
of gas fueling. Unfortunately, the spin measurements through 
these two techniques at this stage seems incompatible to allow 
us draw a firm conclusion. We find that generally a lifetime of 
gas fueling as A7f = 1yr can reproduce the spin distri¬ 
bution reported by the Iron Ka line fitting, whereas a lifetime 
as ALf = 10 5 Mg 5 yr can reproduce these constraint through 
the radiative efficiency. However, the later case seems con¬ 
sistent with the observed episodic lifetime as long as ~ 10 7 yr 
for very luminous AGNs at high redshift, which plausibly har¬ 
bor > 10 9 Mq SMBHs. Since the lifetimes of gas fueling are 
linked to the episodic lifetimes of AGNs, we proposed that 
the episodic lifetimes should be regarded as a new ingredi¬ 


ent for the semi-analytic models of SMBH growth and spin 
evolution. 
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